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Abstract 

We apply some recent developments of Baldoni-DeLoera-Vergne [J on vector partition functions, to 
Kostant and Steinberg formulas, in the case of A r . We therefore get a fast Maple program that computes 
for A r : the multiplicity of the weight p in the representation V(X) of highest weight A; the multiplicity 
c A.ji °f the representation V(v) in V(X) ® V(p). The computation also gives the locally polynomial functions 
C\ and c^ M . 

1 Introduction 

In this short note, we are interested in the two following problems in the case of A r : 
Mult: Computation of the multiplicity of the weight p in the representation V(A) of highest weight A. 

Tens: Computation of the multiplicity of the representation V(y) in the tensor product of representations 
of highest weights A and p,. 

The approach to these problems is through vector partition functions, namely number of integral points 
in lattice polytopes. More precisely, let <E> be a n x N integral matrix with column vectors $1, . . . , <3> at - Fix 
a n-dimensional vector a. The rational convex polytope associated to $ and a is 

P($,a) = I as € R N ; ^ = a, i, > o| . 

We assume that a is in the cone C($) spanned by non-negative linear combinations of the vectors <E>i. We 
also assume that ker(<E>) n R+ = {0}, so that the cone C($) is acute. The vector partition function is then 
by definition 

fc($,a) = |p($,a) nn N \ , 

that is the number of non-negative integral solutions (xi, . . . ,xm) of the equation X^iLi x i®i — a - If $ is the 
matrix of positive roots of A r , then fc($, a) is denoted by k(A^ , a) and it is called Kostant partition function. 
Note that $ is the (r + 1) x (r(r + l)/2) matrix with columns a — ej (l<i<j<r + 1), where e* is the 
canonical basis of R r+ . 

Let £r+i be the set of permutations of (r + 1) elements. This is the Weyl group of A r . Kostant multiplicity 
formula asserts that 

< = ^2(-l) e(w) k[A+, w(\ + p)-(n + p)) , (1) 

w 

where p is half the sum of positive roots. Here, the sum is over the elements w £ E r +i such that w(X + 
p) — (p + p) is in the cone generated by non-negative combinations of positive roots. Moreover e(w) is the 
signature of w. 

Steinberg formula asserts that 

c^, M = Yl (-iy {m)£(w,) k(At,w(\ + p) + w'(p + p)-(v + 2p)) . (2) 

w,w' 

Here, the sum is over couples (tu, w') G S r +i x E r +i such that w{\ + p) + w'(fi + p) — {v + 2p) is in the cone 
C(A+). 

We use results of Baldoni-DeLoera-Vergne [Q and Baldoni-Vergne : 2; on vector partition functions to 
obtain an efficient Maple program. Vector partition function is computed via inverse Laplace formula, 
involving iterated residues of rational functions. 

Recall that LIE program (see [S]) uses Freudenthal and Klymik formulas. The program LIE is designed to 
work for any root system, while our program is designed specially for large parameters in A r . 
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2 Baldoni-DeLoera-Vergne formula 



Consider a r + 1 real dimensional vector space. Let At (the positive root system of A r ) be defined by 

At = {(e l ~e ] );l<i<j <(r + l)}. 
Let E r be the vector space spanned by the elements (e, —ej). Then 

E r — {a £ R r+1 ; a = aiei + ■ ■ ■ + a r e r + a r +ie r +i with a± + a2 + ■ ■ ■ + a r + a,--|-i = 0}. 
The vector space E r is of dimension r and the map 

/ : R r — > £V (3) 

defined by 

a = (ai, a2, ■ • ■ , a r ) 1 — ► ffliei + ■ ■ ■ + a r e r — (oi H h a r )e r+ i 

explicitly provides an isomorphism of E r with the Euclidean space R r . The hyperplane arrangement (setting 
z r +i — 0) generated by At is given by the following set of hyperplanes in C r : 

{zi ■ 1 < i < r} U {(Zi -Zj);l<i<j< r}. 



Let Ra t be the set of rational functions f(zi, z 2 , 



on C r , with poles on the hyperplanes Zi = Zj or 



= 0. For a permutation iw £ S r define the linear form on Ra,, 

IRes™ =0 / = Res Zm(1)=0 Res Ztu(2) = ■ ■ ■ Res*„, (r) =o f(zi, 22, ■ ■ ■ , z r ) 

= Res zl=0 Res Z2 = • ■ ■ Res Zr = f(z w -i w , z w -i (2 ), ■ ■ -,z w -i( r ))- 

In particular for w — id the linear form / 1— » IRes z= o / defined by 

IRes z= o / = Res zl= o Res Z2=0 ■ • • Res Zr = f(zi,z 2 , . . . , z r ) 

is called the iterated residue. 

Let C(A~£) C E r be the cone generated by positive roots. A subset a of A^ is called a basic subset if {a} 
form a vector space basis of E r . The chamber complex is the polyhedral subdivision of the cone C(A r ') which 
is defined as the common refinement of the simplicial cones C(a) running over all possible basic subsets of 
At. The pieces of this subdivision are called chambers. See Q and for the computation of chambers for 

A wall is a hyperplane in E r spanned by (r — 1) vectors of At- A vector v 6 C(At) is regular if it is not 
on a wall. This means that for every strict subset / C {1, . . . , r + 1} we have J^iel Vi ^ ^' 



Let Sp(a) be the set of permutations w £ E r such that: 



if a, 



(i) 



> 



then w(l) < w(2) 
then w(2) < w(3) 



else w(l) > u?(2) 
else w(2) > w(3) 



if a u ,(i) + a w (2) > 

if a w ti) + ■ ■ ■ + a w u) > then w(i) < w(i + 1) else w(i) > w(i + 1) 

if a w n\ + ■ ■ ■ + ffltu(r-i) ^ then w(r — 1) < w(r) else w(r — 1) > w(r) 

An element of Sp(a) will be called a special permutation for a. Remark that if ai > for all i < r, then 
Sp(a) = {id}. Also remark that Sp(a) is a subset of the subgroup E r of the Weyl group E r +i of A r . 
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Given a £ C(At) PI Z r+ , define def(a) — a + e(X^i=i e » ~ re r +i) with e : 
Lemma 2.1 TTie deformed vector def(a) verifies: 

• def(a) is regular. 

• a £ C(A+) i/ and on/y i/def(a) £ C(At). 
Proof: 

• Let 7 a be a strict subset of {1, . . . , r + 1} such that X^ igJ def (o), = 0. 
First, assume that r + 1 ^ 7 a . We can re-index a in order to get I a = {1, . . 
(ai +£) + ■■■ + (0^ + e) = means that the integer ai + ■ ■ ■ + dk equals — 
< 27 < 1, contradiction. 

Now, assume that r + 1 £ I a - We can also assume that I a — {1, . . . , k, r + 1} with k < r. By definition 
(ai +£) + •■• + (a*, + e) + (a r +i — re) equals 0, therefore the integer ai + ■ ■ ■ + at + a r +i is equal to 
But k < r leads to — | < < i, hence k = r. Consequently 7 a = {1, . . . , r + 1}, contradiction. 

• Note that the coordinates ai of a are integers. Now the integer ai + ■ • • + a* is non-negative if and 
only if ai + • • • + Oj + 5- is non-negative, because < ^7 < 1. Hence a £ C(At) is equivalent to 
def(o) £ C(A+). 



, k} with k < r. Thus 
But < k < r implies 
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Now we can state the formula that was implemented: 

Theorem 2.2 (Baldoni-DeLoera-Vergne 1 ) For a g C(A^) n Z r+1 , the Kostant partition functic 
given by: 

k(At,a)= V (-!)"(") IRes^ f Q ± 1^1+ *g 

»6S P (a') V 2: l---^lll<i<j<r(^ 

z J a if a is regular, 
\ def(a) otherwise. 

In particular, if ai > /or 1 < i < r, we have 



k(A~r, a) — Res zl= o Res Z2= o • • ■ Res Zr= o 



3 Deus ex machina 



(1 + Z 1 )° 1 + r - 1 (l + Z 2 ) a2+r - 2 ■■■(! + Z r )° 
Zl - Z r Y[ 1<i<j<r {Zi-Zj) 



This section features a brief description of the algorithms that were implemented with the software Maple. 
This program is available at http://www.math.jussieu.fr/~cochet 

3.1 How to use the program 

The initial data are only vectors: two for computing the multiplicity c^, three for computing the tensor 
product coefficient c\ . 

Our program works with weights represented in the canonical basis of R r+1 , and not fundamental weights 
basis of A r like LIE. The translation between these two approaches is performed via the procedures fundamental 
and fundamentaLinverse. For example fundamental([2,l,-3]) returns [1,4]. 

Therefore computing the multiplicity is done by typing in multiplicity(lambda,mu) where A and fi are 
lists of r + 1 rationals such that X^=i ~ Ei=i M» an d Aj — Aj+i G N, fii — y-i+i G Z. 

For computing the tensor product coefficient c\ , the syntax is tensor_product(lambda,mu,nu) where A, 
H and v are lists of r + 1 rationals such that Yli=i (^» + A 4 *) = v i an( A ^« — •Vfi £ N, fj,i — /Ui+i G N, 

^: — fi+i G N. 

In the examples, we use the vector 8 — rei + (r — l)e2 + • ■ ■ + le r _i — r ( r + 1 ) e r +i- Its decomposition in 
the fundamental weights basis is the r-dimensional vector (1, ...,1,1 + r(r + l)/2). 

3.2 Implementation 

The elements we need to compute are: 

1. The vector a' = def(a) obtained by deforming the initial parameter a. 

2. The residues that appear in theorem 12.21 

3. The two sets of permutations that appear in Kostant and Steinberg formulas (see (0 and 

4. The set of special permutations Sp(def(a)). 

Because of lemma l2~T1 we may use def(a) instead of a and we do this to simplify the procedures. We 
compute the vector def (a) via the straightforward Maple procedure defvector. This takes care of the first 
part. 

Computation of residues is done iteratively. The function F which residue we need to compute is a product 
of a certain number of functions. This allows to take the residues by introducing little by little the part of 
the function F containing the needed variable. See a detailed explanation of this procedure in 

Let ii,i)6 E r . A valid permutation for u and v is a permutation w G £ r +i such that w(u) — v G C(A^). 
We denote by V(u, v) the set of valid permutations for u, v. Hence, Kostant formula for A r rewrites as 

<%= £ (-i) e(w) KA+, w(\+ P )-(jM + p)). 

Given a set of chambers {C w } we v(x,u,) of C(A^), it follows from |7| that is polynomial when wX — fj, £ 
C m , for w G V(\,fi). In particular, the function N i— > C]va,jv m is a polynomial in N of degree less of equal to 

r (r — 1) 
2 ' 

Let us explain our implementation with the symbolic langage Maple of the procedure valid .permutations 
designed to find the set V(u,v). The method is quite simple: we build the permutations iteratively. This 
allows us not examining all permutations and saving much time. Recall that we have to find all permutations 
w's such that u TO m > vi, U^m + Mmp) > t>i + V2, etc. For any sequence x of indices, we denote by u x the 
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Step 1. Let X be the set of all indices i such that Ui > vi. 

Step 2. For each x G X, we find all indices ia; such that u x + Ui m > v x + Vi m . Let X new be the set of such [a;, i x ], 
for all x and i x . Then X <— X new . 

We repeat r times step 2, and obtain the list X of (r + l)-uples representing permutations of 1, . . . , r + 1. 
The second step is treated in the procedure next_index_valid_permutations. The procedure valid_permutations 
contains first step and a for. . . do loop executing r times step 2. 

Remark 3.1 We reduce computing time by using the following three tricks. 

1. We compute once and for all the vector v' — [vi, vi + V2, ■ ■ ■ , vi + • • • + tV+i]- 

2. We build at the same time of X — [[ii, . . . ,i p ], . . . , (other sets of indices)] the set SX of partial sums 
associated to each [ii , i p ] . More precisely SX = [ui 1 + ■ ■ ■ + Ui p , . . . , (other partial sums)] . 

3. We use tables instead of lists. 

Now let us examine the couples of permutations involved in Steinberg formula. Let iti, 112, v € E r . A valid 
couple of permutations for u\ , U2 and v is a couple (wi, W2) £ Sr+i X £r+i such that wi(ui) + 102(112) — v £ 
C(A^). We denote by V(u%, U2, v) the set of valid couples of permutations for ui, U2 and v. Hence Steinberg 
formula rewrites as 

c^= J2 (~lT (wMw,) k(A+,(w(\ + p)+u,'(^ + p)-(u + 2p)). 

(w,w')£V(\+p,n+p,u+2p) 

The procedure computing valid couples of permutations is similar to the former. 

To compute the subset Sp(a) of E r , we use the procedure speciaLpermutations. This procedure is very 
similar to the previous one. We stress that the Maple function combinat[permute] is impractical and does 
not go very far because of memory limitations. 



4 Test of the program 

Let 6 be the r-dimensional vector (1, . . . , 1, 1 + r(r + l)/2) (fundamental weights decomposition in A r ). It 
translates as (r, r — 1, . . . , I, — r(r + l)/2) in the canonical basis of R r+1 . We used this vector to check the 
well-known fact that the multiplicity of the weight in the representation of A r of highest weight NO is given 
by the dimension of the representation of A r -i of highest weight Np, which is (N + l) r ( ? '~ 1 '/ 2 . 

In this test, we compute for various A r (r = 1, . . . , 8): 

• with A = N9 and either fi = (worst case), or fi = [9N/10]9 (intermediate case). 

• c a,m w ith A = fJ, = NO and either v = (worst case), or v — 2[9N/10]8 (intermediate case). 

Tests were made with bi-processor PHI 1, 13GHz. The notation "— " in an array means that we did not 
try the computation (and not that computation failed). 

Recall (see for example [1] and 0) that counting integral points in a lattice polytope is polynomial in the 
size of input if dimension is fixed, and NP-hard if dimension is not fixed. The figures and [5] emphasizes this 
result. 

In figure the letter I stands for intermediate case (fi = [9N/10]9 for and v = 2[9iV/10]# for c|( M ), 
while W stands for worst case (p — for and v — for c\ ^). 





N = 10 1 


N = 10 2 


N = 10 3 


N = 10 4 


N = 10 5 


N = 10 e 


N = 10 7 


N = 10 s 


N = 10 9 
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52.6 s 


53.8s 


57.0s 


58.4s 


59.9s 


62.3 s 


Ca,„> W, A 4 


13.5s 


13.7s 
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Figure 1: Time of computation, when size of input grows 



4 



The computation can also be done with parameters, giving (N + l) r ' r as expected. 
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Figure 2: Multiplicity of in V(N9) when rank increases 
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